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CN . Abstract 

' This paper presents a global optimization approach to quantum mechanics, which describes the most 

| fundamental dynamics of the universe. It suggests that the wave-like behavior of (sub)atomic particles 
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systems can find their ground states corresponding to the global minimum of some energy function 
associated with the system. The classic time-independent Schrodinger equation is shown to be derivable 
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from the global optimization method to support this argument. 



I. Introduction 

Could quantum dynamics follow a global optimization process so that (sub)atomic systems like molecules 
can most often evolve into their ground states? Nature possibly demands such a process so that (sub)atomic 
systems can find their ground states corresponding to the global minimum of some energy function 
associated with the system. If the dynamics of any (sub)atomic system make it evolve into a final 
state which is the global minimum of its associated energy function, then the universe is likely to be 
deterministic. Otherwise, if a local minimum is often reached instead, which can be any one of many 
local minimums of the energy function, then it is likely to be probabilistic. The dynamics which leads 
to deterministic outcomes for the evolution of (sub)atomic systems would make the world more stable. 

Classic optimization methods such as gradient descent are usually based on local optimization to find 
a local optimum instead of the global one (see Chapter 5 in [8] for a comprehensive survey of the 
classic methods for nonlinear programming). The fundamental principle of those methods is based on the 
iterative improvement of solutions which suffers from the local optimum problem. The key characteristic 
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of their operations is to assign each variable a precise value at any given time. If nature deploys a local 
optimization process to define the dynamics of atomic systems where each object (particle) in the system 
has a precise location (or momentum) at any given time, then molecules including proteins will unlikely 
be stable because their final states could be probabilistic and unpredictable. The system can be trapped 
into one local minimum or another, sensitive to initial conditions and perturbations. 

Cooperative optimization is a general optimization method for finding global optimal solutions instead 
of local ones [9], [10], [11], [12], [13]. It does not struggle with local minima and offers us a complete 
departure from conventional local optimization methods. It was originally developed for solving real- 
world, NP-hard optimization problems raised from the areas such as communications [14], computer 
vision [10], [12], and image processing [15]. In real-world applications, cooperative optimization is 
remarkably successful and has often outperformed state-of-the-art algorithms. The key characteristic of 
its operations is to make soft decisions at assigning variable values at any given time; namely, the 
preference of picking each value for a variable is weighted numerically. This paper demonstrates that 
one of the most fundamental equations in physics, the time-independent Schrodinger equation [1], [2], 
[3], can be derived from the dynamic equations of cooperative optimization. Such a connection between 
the quantum dynamic equation and a global optimization algorithm offers us new insight into quantum 
mechanics. 

We postulate that quantum dynamics follows a global optimization process such that (sub)atomic sys- 
tems can most often find their ground states. The wave-like behavior of (sub)atomic particles, besides Max 
Born's probabilistic interpretation [4], could also be the critical characteristic of the global optimization 
process making soft decisions at assigning variable values at any given time so that the global minimum 
instead of a local minimum can be found. 

II. Cooperative Optimization For Finding Global Optima 
Cooperative optimization is a general principle for finding global optima of a multivariate function 
E(x\, X2, ■ ■ ■ , x n ). It utilizes a function in simple forms, such as ip( x ii x 2, ■ ■ ■ , x n ) = J2i^( x i)' an ^ 
iteratively refines the function ^(x) as the lower bound of the multivariate function E(x). At given time 
t, if the lower bound function ^(x,t) has been tightened enough so that its global minimum equals to 
the global optimum of E(x), i.e., 

min^>(x,i) = min£'(x) , 
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then the global minimum of E{x) is found which is the same as the global minimum of &(x,t), 

&TgmmE(x) = argmin ^(x, i) if minf^i) = minE(x) . 
The global minimum of ty(x,t) of the form £V ^(xi,t) can be easily computed as 

x*(^>(x,t)) = argmm^i(xi,t) for i = 1, 2, . . . ,n . 

Xi 

Assume that the multivariate function E[x), often referred to as the energy function in physics, can 
be decomposed as the aggregation of a number of sub-energy functions, 

E{x) = Y,E l {x) . 

i 

Assume further that ^f(x, t) of the form ]TV Vl/(a^, t) is a lower bound function of E(x), ty(x, t) < E(x). 
Let ty(xi, t + 1), for i = 1, 2, . . . , n, be computed as follows 

%( Xi ,t + l)= min \(l-X)E i (x)+Xy^w ij ^ j (x j ,t)\ . (1) 

XjGX^Xt I j J 

Then the new function ^f(x, t + 1) = ^(xj, t + 1) is also a lower bound function of E(x), guaranteed 
by theory [11]. In the above equation, A and w-ij (1 < i,j < n) are coefficients of the linear combination 
of Ei and tyj(xj,t). X satisfies < A < 1 and u^s satisfy Wij > and ^ Wij = 1 for 1 < i,j < n. 

The difference equations (Q define the dynamics of cooperative optimization. The original minimization 
problem min x E(x) has been divided into n sub-problems of minimization (see Q)- Those sub-problems 
can be solved in parallel in implementation. The function ^i(xi,t) is the solution at solving the ith sub- 
problem. The energy function of the ith sub-problem, denoted as Ei{x), is a linear combination of the 
original sub-energy Ei{x) and the solutions from solving sub-problems other than the zth, i.e., 

Ei{x) = (1 - A) Ei(x) + A mj^j(xj,t) . 

j 

The cooperation among solving those sub-problems of minimization is thus achieved by having each sub- 
problem compromising its solution with the solutions from others. Ei(x) is called the modified objective 
function for the sub-problem i. 

The coefficient A is a parameter for controlling the cooperation at solving the sub-problems and is 
called the cooperation strength. A high cooperation strength leads to strong cooperation at solving the 
sub-problems while a lower cooperation strength leads to weak cooperation. The coefficients u>ij control 
the propagation of the sub-problem solutions in the modified objective functions Ei(x) (details in 
[11])- 
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The function ^(xj, t) can be understood as the soft decision at assigning the variable X{ at minimizing 
Ei(x). The most preferable value for variable x% at time t is argmin.^ \E^(xj,t). 

The cooperative optimization theory [11] guarantees that the dynamic system defined by the difference 
equations (0 has a unique equilibrium and converges to it with an exponential rate regardless of initial 
conditions and is insensitive to perturbations. 

Without loss of generality, assume that all energy functions including the sub-energy functions Ei(x) 
are nonnegative functions. Then the cooperative optimization theory tells us that the lower bound function 
^ computed by {TJ can be progressively tightened, 

V(x, t = 0)< *(x, t = l)< *(x, t = 2) < . . . < $(x, t = k) < E(x) , 

when we choose the initial condition as VPi(xi, t = 0) = 0, for i = 1,2, ... ,n. 

With certain settings of the cooperation strength A and the solution propagation coefficients Wij, 
the time-independent Schrodinger equation can be derived in mathematical form from the difference 
equation Q of cooperative optimization. The following two sections offer the detail of the settings and 
the derivation. 

III. Variations of Cooperative Optimization 
If the energy function E(x±,X2, ■ ■ ■ , x n ) is of the following form 

n n n 

E(xi,x 2 ,...,x n ) =y^ei(xj) +y~] eij(xi,Xj), (2) 

i=l i=l j=l,j^i 

Let the decomposition of E{x) be 

Ej = ej(xj) + ejj(xj,Xj), for i = 1,2, . . . ,n (3) 

The difference equations (Q become 

^i(xi,t + 1) = (1 - X)ei(xi) + \wii^i(xi,t) + min((l - X)eij(xi,Xj) + Xw^ j(xj,t)) . (4) 

Choosing the coefficients as u;^ = and Wij = a (a positive constant), for i / j, we have 

^•(xj, t + l) = ei(xi) + ^ min fe^Xj, Xj) + a^'^Xj, t)j , (5) 

where a = Xa and \f^(xj, i) = ^(x^, i)/(l— A). In (|5jl, the parameter a controls the cooperation strength. 
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To possibly improve the convergence of the cooperative optimization defined by the difference equa- 
tions ©, we can offset both sides of the equation by a value, denoted as Z{(t + 1), for i = 1, 2, . . . , n. 
One possible choice for Z{(i + 1) is 

Zi(t+1) =mm^' i (x i ,t + l) . 

Xi 

Let S&ifei, t + 1) = ^S^Xi, t + 1) — Zi{t + 1), the difference equation (J5J) becomes 

^'■(xj,t + 1) = ej(sj) + ^ min (ejj(xj,Xj) +a^(xj,t)) —Zj(t + l) . (6) 

Dividing the both sides of the above equation by a negative value —h (h > 0) followed by taking the 
exponent of the both sides, we have 



(V) 



where 

^(x h t + l) = e-<^ t+1 ^ h . 

The function ipi(xi,t) is the soft decision at assigning the variable x%. It is called the soft assignment 
function in this paper. For each value of X{, the function ipi(xi,t) measures the preference of assigning 
that value at the iteration time t for that variable. A high function value represents a strong preference 
of assigning that value for the variable. 
Using the approximation 

g-max.jfx)//! ^ e -9{x)/h 
x 

which has the property of 

lim e" maX:c 9 ^y h I e -9(p)/n = i 

X 

and choosing the parameter a = 2, the difference equations Q can be rewritten as 

^( X ,,t + 1) = -J_ e -*(«0/S J] [^e-^^)/^,.^,*)! 2 ) , (8) 

where Z»(t + 1) = e*^ 1 )/* 

If t/jj(xj,t) 6 C (the complex domain), the difference equation © still works for global optimization. 
In this case, i)| 2 = tpj(xj,t)ipj(xj,t) represents the soft decision at assigning the variable The 

best candidate value for assigning Xj at time t is the one of the highest function value \tpj{xj,t)\ 2 . Any 
variable value other than the best one may also have a positive function value \ipj(xj,t)\ 2 representing 
the degree of preference of assigning that value to the variable Xj. If only one value of Xj, say x*, has a 
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positive function value \tpj(x*,t)\ 2 > while all other variable values have \ipj(xj,t)\ 2 = (xj ^ xlj), 
then the decision at assigning Xj becomes a hard one instead of a soft one. 

As discussed above, given any i and t + 1, Zi(t + 1) can be any value that can possibly improve the 
convergence. One way to do that is to choose a value for Zi(t), which is equivalent to choose a value 
for Zi{t), such that ^T, x . \^i{ x ii *)| 2 * s bound to a constant, say 1. To be more specific, we choose Zi(t) 
as 

With such a choice, 

^|^(xi,t)| 2 = l . 

Xi 

Hence, this choice makes \ij)i\xi,t)\, for i = 1,2, ... , n, have the property of a probability density 
function. Zj(t + 1) is thus called the normalization factor. (The specific choice of Zi(t + 1) makes no 
difference at the optimization power of the cooperative optimization.) 

If all variables XjS are in a continuous domain, the difference equations © become 



(9) 
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IV. A Continuous Time Version of the Cooperative Optimization algorithm 

Let the soft assignment function of Xi at time t be ipi(xi,t). Let At be an infinitesimal positive 
value and the soft assignment function at t + At be ipi{xi,t + At). The difference equations © of the 
cooperative optimization algorithm in a continuous time version become 

A(xi,t + At) = _L_-^( a . il t)e-( A */ R ) e '( a '«) n / dxj e-^/^^l^^.t)! 2 . (10) 
From dTOb we have 

Yim^ 4>i(xi, t + At) = 4>i(xi,t) . 

A dynamic system described by dTOb is a dissipative system, not a conservative system. It will evolve 
toward its equilibriums over time. It will be shown in the following that the dynamic equation dTOb at its 
equilibrium is, in fact, the time-independent Schrodinger equation. 

Assume that the initial condition is ipi( x 'h *o) = ~ a i)> for i = 1, 2, . . . , n, where 8{xi — ai) is the 
delta function defined as 

S( Xi - ai )=0, when a, and / Sfr - a,) dx, = 1 . 
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and ai is not the optimal value, <n ^ x*. Then (^(xj, to)) is a undesired stationary state of the difference 
equation system (fTOt . To improve the performance of the system, we spread the soft assignment function 
by a smoothing kernel K{x), i.e., 

j K(u - Xi)ipi(u) du =>■ ipi(xi) . 
If Xi is in the one dimensional space R, we can choose the following Gaussian function as K[x) 

K{x) = J_ - e -* 2 / 2 °> At . (11) 
With the soft assignment function spreading, the dynamic equation dTOb becomes 

*><-«- t+ "> - ^^/ d "TO e " ( """' ,V2 '' A '* (M)e " <a,/ * )e ' w 

x H / dxj e-WV^'^jix^t)] 2 , (12) 

Expanding the right side of the above equation into a Taylor series with respect to At and let At — > 0, 
we have the new differential equations for the cooperative optimization algorithm in a continuous time 
version, 

dibi(x,t) af d 2 ib(xi,t) TT , .1 . , . ,, N , . . „ , . 

dt = ~2 dx 2 ~ ^)^(^) + £ i( t )M*i>t), for each (13) 

where 

d Zi(t)/d t 



Vi{xj) = ej(xj) + ^ J dxj eij(x i ,Xj)\'4>j(xj,t)\ 2 and Si(t) 



Let the operator V? be defined as 



X72 ,f + \ d 2 ii(xi,t) 



dx 



and Hi be an operator on ip(xi,t) defined as 



fli = -^V? + Vi(xi) . (14) 

Then the equation (fT3l can be rewritten as 

^i^lH = -lHiil>i(xi,t) + £i(t)i(>i(xi,t) . (15) 
at n 

When the differential equations H51 evolve into a stationary state (equilibrium), the differential equa- 
tions ( fTSI l become 

Eiipi(xi,t) = Hii/ji(xi,t), for i = 1,2, . . . ,n , (16) 
where Ei, E\ = fej, is a scalar. 
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For a physical system consisting of n particles, let X{ be the position of particle i, 1 < i < n, in the 
one dimensional space M. Let 

of = h/nii, 

where rrii is the mass of particle i. Then equations dT6t become 

EiiPi(xi,t) = (-^T V ? + ^(^)) Mxi,t), for i = 1,2,..., n . (17) 

They are the condition for the physical system being in a stationary state when its dynamics is defined 
by the cooperative optimization. Equation H71 is also the time-independent Schrtidinger equation. (It 
is straightforward to generalize this derivation to three dimensions, but it does not yield any deeper 
understanding.) From the equation (fT7l we can see that the soft assignment function ipi(xi,t) is the 
wavefunction in the Schrodinger equation. 

If we choose more advanced forms for ^(x) where two (or more) variables, say Xi and Xj, are clustered, 
together with their soft assignment functions ^i(xi) and S&j(xj), i.e., ^(xj), ^j(xj) =^ ^(xj, Xj), then 
better lower bound functions in terms of tightness can be found and there is a higher chance for the 
cooperative optimization algorithm to find a global minimum [10]. However, more variable clustering 
leads to higher computational cost as a tradeoff. The quantum entanglement phenomena could possibly 
be understood as variable clustering in cooperative optimization to increase the chance of finding the 
global optimum. 

V. Conclusions and Future Research 

This paper derived the time-independent Schrodinger equation being as the dynamic equation of coop- 
erative optimization at its equilibrium in a continuous-time version for continuous variables. Cooperative 
optimization is a general method for finding global optima with the key characteristic of soft decision 
making in assigning variables. The soft assignment functions used by cooperative optimization can be 
interpreted as the wavefunctions in the Schrodinger equation. It could be a critical feature for (sub)atomic 
systems to have deterministic final states, namely their ground states. 

In classic quantum mechanics, wavefunction spreading is understood to have the purpose of preventing 
electrons from collapsing into the nucleus caused by the Coulomb force. In the global optimization 
approach, it also serves the purpose of improving convergence property of quantum dynamics at finding 
global optima. 

For a closed system, we can use the time dependent Schrodinger equation to describe its evolution 
(nonrelativistic). Usually, a closed system is a conservative system. However, it is hard to find a closed 
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system in nature because the (sub)atomic particles in the system will inevitably interact with their 
environment and the system becomes dissipative. Could global optimization be an emerging property 
of such a dissipative process so that (sub)atomic systems can most often find their ground states? Could 
the wave-like property of (sub)atomic particles play a key role at defining such a global optimization 
process? Could the dynamic equation offered by cooperative optimization be suitable for describing the 
global optimization process at a proper level? These are interesting questions worth pursuing with future 
research. 
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